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DISTINGUISHING NON-RESONANT 
FOUR-WAVE-MIXING NOISE IN COHERENT 
STOKES AND ANTI-STOKES RAMAN 
SCATTERING 

5 

FEDERALLY SPONSORED RESEARCH OR 
DEVELOPMENT 

The subject matter of this application may have been 
funded in part under Contract No. NAS2-02057 awarded by to 
the National Aeronautics and Space Administration (NASA). 
The U.S. Government may have rights in this invention. 

BACKGROUND 

15 

Molecules frequently have molecular resonance frequen- 
cies that are due to the electromagnetic attractions of atoms in 
the molecule. These frequencies are those of molecular vibra- 
tions, molecular rotational motions, the excitation of elec- 
trons to higher energy states, and occasionally finer structures 20 
such as hyperfine interactions and optical-magnetic proper- 
ties. These properties are present without the introduction of 
any external contrast molecule. These frequencies are usually 
in the mid-infrared, corresponding to photons of 1.5-50 
microns of wavelength. Unfortunately, they cannot be 25 
directly excited by electromagnetic radiation of the same 
frequency because when they are in tissue, the surrounding 
water absorbs almost all of these frequencies. The range of 
wavelengths that the tissue is relatively transparent to is 0.6- 
1 .5 microns. Therefore multiphoton nonlinear processes need 30 
to be employed to probe these resonances. The photons to 
stimulate and record the processes are typically in a region 
where the tissue is not absorbing, so that they can reach the 
tissue feature and be measured from the feature. 

Raman spectroscopy, first discovered in 1928, uses 35 
molecular resonance features of frequency Act) to split a pho- 
ton of frequency 00 into another photon of frequency co-Ao) 
and a resonance excitation of frequency Aco. The presence of 
photons at frequency co-Aco identifies the concentration of the 
resonance feature. This process is in practice very weak and 40 
requires large amounts of power to produce any detectable 
amount of photons. This weakness is due to the fact that the 
probability of a Raman excitation process to occur is propor- 
tional to the number of photons at frequency co-Aco already 
present, of which there are typically few or none. Since pho- 45 
tons that would be emitted by Raman excitation at frequency 
co-Aco are indistinguishable from the incoming radiation that 
stimulates them, this is not a viable technique for achieving 
molecular sensitivity. 

Coherent Anti -Stokes Raman Scattering (CARS) is 50 
another nonlinear spectroscopy technique that unlike conven- 
tional Raman spectroscopy, allows all of the photons neces- 
sary to stimulate the process to be introduced into the tissue 
by the illuminating source. This enables the probability of a 
CARS interaction to be increased to a (theoretically arbi- 55 
trarily) high level so that a sufficient number of photons can 
be produced as to enable detection within a reasonable time 
period. It is essentially two stimulated Raman scattering pro- 
cesses in parallel. Two photons, a “pump” of frequency 00 x 
and a “Stokes” of frequency oo 2 illuminate the tissue. They 60 
must be separated in frequency by co 1 -cd 2 =Aco, which is the 
frequency of the molecular resonance. When molecules of the 
target molecular species are present, the resonance will be 
excited, and the pump photon will be converted to the same 
frequency as the Stokes photon. This is the first stimulated 65 
Raman scattering process. Another photon may arrive at fre- 
quency oo 3 that will stimulate the emission of the excitation 
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from the resonance, so that the energy of the photon of fre- 
quency co 3 and the excitation are converted to a new photon of 
frequency 004 = 003 + Aco, called the “anti-Stokes” photon. The 
presence of this photon of frequency oo 4 indicates that a 
CARS process has taken place and indeed a molecule with the 
resonance feature is present. Often the “pump” beam is used 
as the photons of frequency oo 3 , so that oo 3 =oo 1 and oo 4 = 2 oo 1 - 
oo 2 . Since the photon of oo 4 is not the same frequency as one of 
the illuminating photons, and is typically within the transpar- 
ency range of the tissue, it is easily discriminated from the 
incoming radiation. FIG. 1 A shows an energy-level diagram 
for CARS, and FIG. 1 B shows an energy-level diagram for 
Coherent Stokes Raman Scattering (CSRS). 

CARS microscopy uses the CARS process to look for the 
presence of a molecular species, but does not require any 
foreign substances to be introduced into the tissue. It scans the 
illumination point-by-point through the tissue and measures 
the number of generated anti-Stokes photons. When a three- 
dimensional mesh of points has been scanned, a complete 
three-dimensional picture of molecules of that resonance can 
be shown. Since CARS is a nonlinear process (and therefore 
is intensity sensitive), efficient conversion only occurs at the 
focus of the illumination, which can be made very tight (typi- 
cally a half micron in both the axial and lateral directions). 
Therefore the resolution can be made many orders of magni- 
tude better than MRI, which is the probably the largest com- 
petition for clinical use for similar purposes. Unfortunately, 
the penetration is usually rather low (limited to about 500 
microns). A further shortcoming is that CARS microscopy 
measures the total number of anti-Stokes photons, or power, 
from the sample. However, the optical field contains temporal 
structure in the phase that is averaged out by power detection 
because photodetector response time is orders of magnitude 
slower than the oscillations of the optical field. The time scale 
on which the optical pulse varies (which is typically picosec- 
onds or femtosecond time scales) is far too fast for photon 
detection equipment or electronics to detect (the fastest of 
which may detect 25 ps time scales). 

Optical coherence tomography (OCT) is an emerging 
high-resolution medical and biological imaging technology. 
OCT is analogous to ultrasound B-mode imaging except 
reflections of low-coherence light are detected rather than 
sound. OCT detects changes in the backscattered amplitude 
and phase of light. 

Nonlinear interferometric vibrational imaging (NIVI) is a 
method used to measure the three-dimensional distribution of 
molecular species in various samples (biological or other- 
wise) [1]. Its basic operation is to stimulate the excitation of 
molecular bonds with particular resonance frequencies, and 
then use these excitations to produce radiation distinct from 
the excitation that can be measured. Unlike previous methods 
that use CARS in microscopy to probe for the presence of 
molecular species, NIVI utilizes a heterodyne approach 
where a reference signal is separately generated and inter- 
feromerically compared to the signal received from the 
sample, allowing the signal to be formed into an image in the 
same way as OCT. In this way, additional information can be 
inferred from the emitted radiation such as the distance to the 
sample and phase information that yields additional structure 
of the molecular bonds. It also has other advantages in sen- 
sitivity and the ability to screen out background radiation that 
is not produced by the sample. It also can allow more flex- 
ibility in the choice of laser illumination source, because the 
coherent detection process does not rely on photon frequency 
alone to discriminate emitted radiation. 

There are compelling reasons to use broadband sources to 
excite CARS. An ultrafast pulse can be shaped into a longer 
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picosecond pulse that can excite CARS more efficiently. In 
addition, it can be used to excite many resonances simulta- 
neously [2]. Unfortunately, unlike narrowband pulses, the 
anti-Stokes radiation produced will not be narrowband. If 
many resonances are excited simultaneously then the anti- 
Stokes radiation they produce have overlapping spectra. 
Because noninterferometric detection can only measure the 
spectrum of the anti- Stokes radiation, the contributions of 
each resonance to the anti-Stokes radiation will be insepa- 
rable. Interferometric detection allows the demodulation of 
the anti-Stokes field so the Raman spectrum can be inferred 
when exciting multiple simultaneous resonances. 

In addition, broadband sources should allow for more 
Raman-frequency agile imaging instruments. When utilizing 
narrowband pump and Stokes pulses, the frequency differ- 
ence between them must be tuned to the Raman frequency of 
interest. Often retuning lasers or amplifiers is difficult to make 
reliable and automatic. Pulse shaping, however, is achieved 
by movable gratings, prisms, or mirrors to adjust dispersion 
and delay, or by acousto-optic or liquid-crystal Fourier plane 
pulse shapers without moving parts [3]. Because pulse 
shapers do not typically involve feedback, oscillation, or 
overly sensitive alignment, pulse shapes can be changed 
much more easily. In addition, a computer can control these 
mechanisms automatically, so that changing the pulse shape 
should be much easier than retuning a laser source. 

Two important features of the Raman spectrum in practice 
that need to be mutually distinguished are the resonant and 
nonresonant components. The resonant components are spe- 
cific to features of a molecule, which include vibrational- 
frequencies, rotational frequencies, and electronic reso- 
nances. Nonresonant features are not specific to a particular 
molecule, and are weakly dependent on frequency. 

Distinguishing CARS resonant from nonresonant four- 
wave-mixing has become problematic with the advent of 
ultrafast laser sources. These sources produce pulses on the 
order of 5-200 fs, much shorter than the lifetime of the reso- 
nance, which is 1/T„. If transform- limited ultrafast pulses are 
used, only a small polarization P ( 3 ) (Q) can be produced in the 
molecule, while the nonresonant components are enhanced. 
Since the lifetime of the resonance is typically 1-100 ps, 
transform-limited ultrafast pulses excite resonant transitions 
inefficiently and nonresonant transitions efficiently. Before 
the advent of near-infrared solid-state femtosecond lasers, 
especially the Ti 3 + :A1 2 0 3 laser, most laser sources produced 
picosecond-length pulses that favored the generation of 
CARS. Many current CARS instruments utilize narrowband 
pump and Stokes pulses of picosecond length for this reason. 

SUMMARY 

In a first aspect, the present invention is a method of exam- 
ining a sample, comprising exposing of the sample to a pump 
pulse of electromagnetic radiation for a first period of time, 
exposing of the sample to a stimulant pulse of electromag- 
netic radiation for a second period of time which overlaps in 
time with at least a portion of the first exposing, to produce a 
signal pulse of electromagnetic radiation for a third period of 
time, and interfering the signal pulse with a reference pulse of 
electromagnetic radiation, to determine which portions of the 
signal pulse were produced during the exposing of the sample 
to the stimulant pulse. The first and third periods of time are 
each greater than the second period of time. 

In a second aspect, the present invention is a method of 
producing an image, comprising collecting a set of data. The 
data comprises a plurality of digital data produced by the 
above method. 
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In a third aspect, the present invention is a method of 
examining a sample by CARS or CSRS, including exposing a 
sample to laser light to produce a signal, and producing a data 
set from the signal. The improvement comprises interfering 
5 the signal with a reference pulse, and excluding at least a 
portion of the signal containing non-resonant electromag- 
netic radiation in producing the data set. 

In a fourth aspect, the present invention is a method of 
examining a sample, comprising a step for producing a CARS 
10 or CSRS signal pulse of electromagnetic radiation from a 
sample, and a step for determining which portions of the 
signal pulse contain electromagnetic radiation produced by 
four-wave mixing. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1. Coherent Anti-Stokes Raman Scattering and 
Coherent Stokes Raman Scattering energy-level diagrams. 

FIG. 2. Basic block diagram of nonlinear interferometric 
20 vibrational imaging (NIVI). 

FIG. 3. Example laser configurations that produce two 
pulses of frequencies oq and oo 2 overlapped. 

FIG. 4. Methods of shaping a broadband pulse into a pulse 
with beat frequency Aoo. 

25 FIG. 5. Reference Signal Generator Implementation 

FIG. 6 . Configurations for full field CARS. 

FIG. 7. Translated serial-point scanning configurations 

FIG. 8 . Beam-steered serial -port scanning configurations 

FIG. 9. Temporal-ranging based NIVI 
30 FIG. 10. Full field cross-correlator demodulator 

FIG. 11. Temporal cross -correlator for a serial -point scan- 
ning microscope. 

FIG. 12. Two configurations that utilize linear photodetec- 
tor arrays to measure multiple samples of the cross-correla- 
35 tion simultaneously. 

FIG. 13. Diagram of an implementation of the recorder. 

FIG. 14. Electrical cross-correlation signal quadrature 
amplitude demodulator. 

FIG. 15. “Biological window” in tissue where absorption 
40 of near-infrared wavelengths is at a minimum and light can 
penetrate deep into highly-scattering tissue. 

FIG. 16. Flowchart of preconditioned conjugate gradient 
algorithm applied to inversion of Raman spectra. 

FIG. 17. Four- wave-mixing (FWM) and CARS signals 
45 produced by a long pump pulse and a short Stokes pulse 
combination. 

FIG. 18. Experimentally measured cross-correlations 
between a short reference pulse and the anti-Stokes radiation 
of an acetone sample as the pump and Stokes frequency 
50 difference are tuned. Away from the resonance at 2925 cm -1 , 
only the nonresonant component exists, while the “tail” is 
maximized when tuned near resonance. 

FIG. 19. Schematic of experimental apparatus used to 
acquire cross-correlation interferograms. 

55 FIG. 20. Setup to generate a CARS excitation signal and 
reference signal to reject nonresonant four-wave-mixing, 
using a delayed probe pulse so that a reference pulse can be 
used as an interferometric gate to reject the early-arriving 
nonresonant four-wave-mixing. 

60 FIG. 21. Setup to generate a reference pulse and over- 
lapped dispersed pump and Stokes pulses for measurement of 

FIG. 22. Setup to measure Raman spectrum in a particular 
frequency range using an ultrabroadband source. 

65 FIG. 23. Graph of instantaneous frequency of two over- 
lapped pulses, where the pulse frequency difference sweeps 
from Q h to Q l . 
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FIG. 24. Time/frequency diagram for pulse combination 
using linearly chirped and delayed pulses. 

FIG. 25. Simulation of measurement and reconstruction of 
Raman spectrum, (a) is the magnitude of the illumination 
pulse as a function of time, (b) is the spectrum of the illumi- 5 
nation pulse, (c) is the spectrum of the generated anti-Stokes 
light, (d) is the spectrum of the beat frequencies of the illu- 
mination pulse, (e) is the magnitude of the Raman spectrum 
of the hypothetical molecule, (f) is the recovered hypothetical 
spectrum. 10 

FIG. 26. Simulation of measurement and reconstruction of 
Raman spectrum, (a) is the spectrum of the beat frequencies 
of the illumination pulse, (e) is the magnitude of the Raman 
spectrum of the hypothetical molecule which includes both 
resonances. (1) is the recovered hypothetical spectrum. 15 

DETAILED DESCRIPTION 

The present invention makes use of the discovery that 
interferometry may be used to differentiate between resonant 20 
CARS and non-resonant four- wave -mixing processes. Inclu- 
sion of a reference signal allows for interferometry, and there- 
fore this important source of noise can be eliminated. This 
may be particularly advantageous when combined with NIVI, 
which already includes a reference signal. 25 

A possible pump/Stokes combination for measuring spec- 
tra using CARS, would include a narrow band pump of fre- 
quency cOj used with a broadband short Stokes pulse of fre- 
quencies oo 2 -Aoo/ 2 to co 2 +Aco/ 2. The Stokes pulse arrives on 
the leading edge of the pump pulse, and excites a broad band 30 
of Raman frequencies (u 2 -co 1 -Aco/2<Q<co 2 -tt) 1 +AeL)/2. 
Because of the narrowness of the pump, a Raman frequency 
of Q will be mapped to an anti-Stokes frequency of oq+Q. A 
spectrometer can then directly sample the anti-Stokes spec- 
trum, which maps onto the Raman spectrum. 35 

While this approach has the appeal of simplicity, it has a 
number of drawbacks that make it difficult to implement in 
practice. The chief problem is generating synchronized nar- 
rowband pump and broadband Stokes pulses. The bandwidth 
of the pump must be equal to or less than the bandwidth of the 40 
Raman spectral lines probed (on the order of 1-10 cm -1 ) 
while the bandwidth of the Stokes must exceed the total 
bandwidth probed (usually greater than 100 cm -1 ). Using 
synchronized picosecond and femtosecond lasers is possible, 
but difficult and cumbersome. Otherwise, both pulses will 45 
usually be derived from a single oscillator, usually a femto- 
second oscillator, with an optical parametric amplifier or 
oscillator used to convert part of the pump pulse to the Stokes 
pulse. A broadband pump pulse must be used to derive a 
broadband Stokes pulse. To decrease the bandwidth of the 50 
pump for this technique, it can be filtered to a narrow band- 
width. However, because the anti-Stokes output power is 
proportional to the square of the pump power the decrease in 
anti-Stokes power will be severe. 

If one attempts to start with a narrowband picosecond 55 
oscillator, and derive the pump and Stokes from that, the 
Stokes pulse will need to be broadened using a process such 
as self-phase-modulation. Unfortunately, these nonlinear 
broadening processes require very high peak power pulses, 
which for relatively long picosecond pulses would require 60 
very high pulse energy indeed. Also, the bandwidth created in 
this way tends to be difficult to recompress to a pulse shorter 
than the original pulse length. 

In addition, the Stokes pulse will need to be chirped to 
avoid overly increasing the nonresonant four-wave-mixing, 65 
which will require further lengthening of the pump pulse in 
time to ensure sufficient Raman frequency resolution. This 
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will further decrease the anti-Stokes output. While sampling 
of the Raman spectrum is possible using incoherent detection, 
it presents difficulties that are much more easily solved using 
coherent interferometric detection. 

A typical Raman spectrum % (3 \Q) can be decomposed into 
a sum of resonance and nonresonant components: 

v 3 i (fl)=^(n)+y / (11) 


The nonresonant component is a slowly chang- 

ing function of frequency Q, and is usually approximated by 
a constants Xnr^ 3 \ < ^) = Xnr which is a real number. The reso- 
nant components x„ ( 3 ) (Q) are sharply peaked at the resonance 
frequency, and are typically described by a homogeneously 
broadened Lorentzian spectrum: 


v/OUm - v ' 

Xn K } £l 2 -2iY n £l-{£ll-Yl) 

where Q=^Q„ 2 - 2 T „ 2 is the center frequency of the reso- 
nance, and T„ is the linewidth. A resonance has a Raman 
magnitude that is sharply peaked around the center frequency. 
Interferometric detection can distinguish the flat, insensi- 
tively dependent on Raman frequency, nonresonant spectra 
and the sharply peaked resonant CARS spectra. In addition, 
the imaginary part of % ( 3 ) (Q) is determined by only the reso- 
nant Raman component, while the real part is determined by 
both resonant and nonresonant components. Another of the 
distinguishing features between a resonant and nonresonant 
feature is the phase reversal of jt that occurs through the center 
frequency. Noninterferometric instruments cannot distin- 
guish the phase, while the interferometric method measures 
the complex 5 C„ ( 3 ) (£ 2 ). As well as being a signature of CARS, 
this phase reversal may be used to separate several resonances 
close together in frequency. 

To see how interferometry can detect the difference 
between resonant and nonresonant processes, consider the 
shape of the anti-Stokes pulses produced by the long pump 
and short Stokes pulses previously mentioned. This combi- 
nation is illustrated in FIG. 17. When the pump and Stokes 
pulses overlap, the molecule will be excited by stimulated 
Raman scattering (SRS). This excitation will remain after the 
Stokes pulse passes. At the moment of the overlap, nonreso- 
nant four- wave-mixing processes can also be excited. How- 
ever, because there is no persistent state associated with non- 
resonant processes, the emission of four-wave-mixing will 
end quickly after the Stokes pulse passes. With a Raman- 
active resonance, however, the pump can produce anti- Stokes 
radiation via SRS even after the Stokes has passed, because 
the resonance persists. Therefore, the nonresonant compo- 
nent can be discarded by rejecting any anti-Stokes radiation 
that occurs coincident with the Stokes pulse. 

The benefit of interferometry is that the time of arrival of 
the anti-Stokes radiation can be found very precisely by 
cross -correlating a broad bandwidth reference pulse with the 
anti- Stokes radiation. Incoherent detection can detect the 
interference between the resonant and nonresonant compo- 
nents of the anti-Stokes radiation, but does not directly detect 
the time of arrival of the anti-Stokes light. FIG. 18 shows the 
actual measured cross-correlations between a short reference 
pulse and the anti-Stokes radiation for a cuvette of acetone 
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with a resonance at 2925 cm -1 as the frequency difference 
between the pump and Stokes is tuned. 

The apparatus to acquire the interferogram is documented 
in FIG. 19. A regenerative amplifier produces a pulse at 808 
nm and 30 nni bandwidth that is used as a seed pulse for a 
second-harmonic -generation optical parametric amplifier, 
and also as the pump for the CARS sample. The pump is 
lengthened to approximately 200 fs by passing it through a 
dispersive Dove prism made of BK7 glass of 105 mm length. 
The idler of the optical parametric amplifier produces a 
Stokes pulse at 1056 nm and 70 fs in length, which is com- 
bined with the pump at a dichroic mirror. The pump is delayed 
so the Stokes arrives on the leading edge of the pump pulse. 
The signal of the optical parametric amplifier is at 653 nm and 
serves as the reference pulse because it matches the frequency 
of the anti- Stokes produced in the CARS sample when the 
pump and Stokes are focused into it. The anti-Stokes and 
reference are cross -correlated by delaying them relative to 
each other and detecting their interference power on a silicon 
photodetector. 

Note that to detect the difference between resonant and 
nonresonant processes did not require that the pump be 
lengthened to the lifetime of the resonance of the acetone. 
This is because the nonresonant component ends abruptly at 
the end of the Stokes pulse. Unlike the incoherent detection 
case, one can actually reject the nonresonant component 
based on its time of arrival, and not just depend on minimizing 
the amount of nonresonant four-wave-mixing that is gener- 
ated. 

Instead of dispersing the pump pulse to create a portion of 
the pump that can act as a probe to produce SRS in the sample, 
the pump pulse can be split into two and one part delayed to 
act as the probe. This can be achieved with the setup in FIG. 
20. A regenerative amplifier or a mode-locked laser oscillator 
generates short pulses which are split by a beamsplitter into 
two pulses. One pulse is to create the pump/probe pulses. This 
pulse is further split in two, and one copy is delayed with 
respect to the other. The earlier pulse is the pump and the later 
pulse the probe. The other pulse not used to create the pump/ 
probe pulses is used in an optical parametric amplifier or an 
optical parametric oscillator to generate nearly transform- 
limited idler and signal pulses. The idler is used as the Stokes 
pulse, and the signal is used as the reference pulse. The idler 
is delayed so that it overlapped with the pump pulse in time 
and space with a dichroic beamsplitter. The pump/Stokes/ 
probe combination is sent to a microscope to excite the 
sample. The probe will create stimulated Raman scattering, 
which is collected by the microscope. The signal is used as the 
reference and is delayed to arrive at the interferometric 
demodulator at the same time as the SRS created by the probe. 
The magnitude of the interference signal indicates the pres- 
ence of the target species in the sample. 

In some cases, the optical parametric amplifier or oscillator 
may not provide a signal at a suitable frequency to use as a 
reference pulse. If a second-harmonic-generation optical 
parametric amplifier (SHG OPA) is used, then the sum of the 
idler and signal frequencies will equal two times the seed 
frequency, and so the signal will be of the correct frequency. 
If a standard optical parametric amplifier or oscillator is used, 
then the signal and idler frequencies add up to the seed fre- 
quency. The signal pulse can be used as the Stokes pulse, but 
the idler frequency is too low to be useful. In this case, one can 
produce a reference by combining nearly transform-limited 
pump and Stokes pulses withfour-wave-mixing. Effectively a 
SHG OPA/OPO implements a four-wave-mixing process 
with two three-wave-mixing processes, but the same can be 
done directly in a nonresonant four-wave-mixing medium 
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such as quartz, sapphire, CaF 2 , LiF, water, or any medium 
with a sufficiently strong four-wave-mixing cross section, 
transparency to the pump/Stokes/anti-Stokes wavelengths, 
and no nearby resonant Raman frequencies. 

5 We now present a method using a pulse sequence like that 
of FIG. 17 to recover the Raman spectrum x (3) (Q) in a par- 
ticular frequency range. The method uses a pump that is 
broadband but dispersed to be lengthened in time, while using 
a short, nearly transform-limited Stokes pulse. These two 
10 pulses will be overlapped so that the Stokes pulse is on the 
leading edge of the pump pulse. FIG. 21 shows a block 
diagram of apparatus that can produce this pulse combina- 
tion: 

The idea of the setup of FIG. 21 is to stimulate a broad 
15 bandwidth of Raman transitions using a short Stokes pulse 
overlapped with the leading edge of a pump pulse. After the 
Stokes pulse is over, the pump pulse continues to produce 
stimulated Raman scattering, which is demodulated with the 
reference pulse. The portion of the cross-correlation that 
20 overlaps the Stokes pulse arrival can be discarded to remove 
the nonresonant four- wave-mixing component of the anti- 
Stokes signal. Because the field used to produce the resonant 
cars is only due to the pump pulse, the cross -correlation 
amplitude and the pump pulse can be used to estimate the 
25 Raman spectrum x (3) (£2) (for example, using an inverse of Eq. 
6 in the Example section, below). 

To elaborate on the setup, a regenerative amplifier or oscil- 
lator generates a train of ultrashort pulses, typically several 
hundred cm -1 in bandwidth. These pulses are split into two by 
30 a beam splitter. One of the pulses is used in a dispersive pulse 
shaper, which imparts a non-constant group delay to the 
pulse. Such a dispersive pulse shaper can be implemented by 
dispersive materials such as optical glasses or other transpar- 
ent materials, prism or grating type pulse shapers, or optical 
35 fibers or waveguides. The pulse will usually be dispersed to a 
length on the order of the lifetime of the resonances under 
study. The other pulse will be sent to an optical parametric 
amplifier or oscillator, where it will generate a broadband 
signal and idler pulse. The idler will act as the Stokes and will 
40 be combined with a dichroic beamsplitter with the dispersed 
pulse at its leading edge. This combination will illuminate the 
sample through the microscope, and will have anti-Stokes 
radiation collected by the microscope. The signal will act as a 
reference pulse and will be cross-correlated with the anti- 
45 Stokes light generated by the sample. For best performance, 
the pump/Stokes signal should be filtered out of the anti- 
Stokes signal before cross-correlation. 

In addition, there is an optional section of FIG. 21 that 
helps enhance the amount of excitation generated in the 
50 sample. Instead of sending all of the pump light into the 
dispersive pulse shaper, it can be further split in two, where a 
portion will enter the pulse shaper and become a probe pulse. 
The other portion will be a pump pulse, and is delayed but not 
dispersed and is recombined with the dispersed pulse. The 
55 undispersed pump pulse will be delayed to overlap the Stokes 
pulse on the leading edge of the dispersed probe pulse. By 
increasing the intensity of the pump at the instant the Stokes 
arrives, the magnitude of the excitation of the resonances can 
be enhanced, allowing a larger stimulated Raman signal to be 
60 collected from the sample. 

A theory of CARS interferometry and inference, including 
Eq. 1-10, is presented in the EXAMPLES section, below. 

To understand in better detail how the present invention can 
recover the Raman spectrum, consider a dispersed pump 
65 pulse and a Stokes pulse with an electric field 

£ p ((o)=^((o) ex p(/ < j)((o)) for (o 0 <co<co 0 +AcOp (13) 
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E p (oj)=0 otherwise 

E s (co)=A s for o) 0 -Q 0 -Aco iS ./2<oo<o) 0 -Q 0 +A(O jS ./2 


(14) 


E s (oo )=0 otherwise 

The frequency Q 0 corresponds to the resonant Raman fre- 
quency of interest, while the bandwidths Ao)p and Aoo s give 
the bandwidths of the pump and Stokes signals respectively. 
The amplitudes A / ,(oo) and correspond to the real, positive 

amplitudes of the frequency components of the pump and 
Stokes signals, respectively. To find these pulses in the time 
domain, we use the stationary phase approximation (which 
should apply well to dispersed pulses). We define d^/doo^t* 
(00), co'(t)=t ,-1 (o)) (the inverse function of f (00)), and dO/dt=co’ 
(t). Note that f (00) must be a strictly increasing or decreasing 
function of 00. The time-domain versions of these signals are: 


\d z <$> 


1-1/2 


E p ( t ) - A p {(o' (r))|-^j- 1 exp(i<&(r)) for f (w 0 ) < t < t* (a> 0 + Aw P ) 

E t (oj) = 0 otherwise 

E s (r) = A s exp(i(w 0 - f)o)r)sin(rAw5)/(rAw 5 ) 


pulse, we assume that d(|)/dcol a 


P (3) (Q)*EpE 5 x (3) (^) for Q 0 -A<a s /2<Q<Q 0 +Ao) s /2 
P (3) (Q )=0 otherwise 

The anti- Stokes generated from this polarization is: 


10 


15 


(15) 


20 


(16) 


To place the Stokes pulse on the leading edge of the pump 


= 0 . At time t= 0 , the Raman 


frequencies Q 0 -Aco 6 ./ 2 <Q<Q 0 + ^ co s /2 are excited. Because 
the Stokes is nearly transform-limited, the phases of Eq. 1 
cancel and the polarization is proportional to the Raman 
spectrum: 


(17) 


35 


10 


in the time domain. The P (3) (Q) and therefore / w (3) (Q) can be 
recovered from P (3) (t) by means of a Fourier transform. 

Raman Spectra Retrieval Using Broadband Pulses 

One of the chief benefits of the invention is that it enables 
the use of ultrabroadband lasers with thousands of cm -1 of 
bandwidth to be used to measure either narrow (less than 10 
cm -1 ) or a very wide range of Raman frequencies. Such lasers 
include ultrabroadband Ti-sapphire lasers, Cr:forsterite 
lasers, dye lasers, Yb:YAG, Yb-silica fiber, Er-silica fiber, 
Nd:glass, and femtosecond lasers (Ti-sapphire, dye lasers, 
Cr:forsterite, Yb:YAG, Yb-silica fiber, Er-silica fiber, Nd:g- 
lass) broadened by continuum generation (e.g. using bulk 
materials, such as photonic crystal fibers, high numerical 
aperture fibers, or micro structured optical fibers, as well as 
sapphire and water). The benefits of this are that the same 
laser can be used to produce light that both acts as the pump 
and Stokes frequencies (rather than requiring a separate 
Stokes to be derived), and the entire power spectral band- 
width of the laser can be used to contribute to the excitation 
and stimulated Raman scattering of the sample. When using 
narrowband pulses, the power spectral density of the pump 
and Stokes must be very high to ensure sufficient excitation. 
With ultrabroadband excitation one can use its relatively low 
power spectral density to excite, but it does not produce a 
one-to-one correspondence between anti-Stokes and Raman 
frequencies, so that interferometric detection is needed. 

The essence behind these methods is a reinterpretation of 
Eq. 1 . In the time domain, Eq. 1 becomes: 


30 


= r 
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(20) 


Eq. 20 is the causal convolution of the impulse Raman 
response 


r +Aoj^/2 

T* 3) (il)A p (o) - n)exp(i0(w - £1)) d£l 
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If we define 


r —Aio/2 

F* 3 \£l)exp(-i£lt)dt 

v Atu/2 


then this becomes a product: 
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E A (t) = P- 3 \t)A p (oj'(t))-^ r exp(»(r)) 
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45 
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(19) 


where E^(t) is the complex analytic continuation of the 
Fourier transform of 64(00). This suggests that P (3) (t) can be 
recovered by multiplying E^(t) by the conjugate of the phase 
of the probe field 


d 2 <s> 

exp (-«$(/)) -jp 


60 


65 


V 3) (0 = 


40 


J y 3) (ri)exp (-i£lt)dt 


with the instantaneous intensity of the incoming pulse. If 
there is a modulation of the intensity at the same frequency as 
a Raman resonance, the Raman resonance will become 
excited. By using a pulse that contains a modulation of the 
intensity of the excitation wave that spans a range of frequen- 
cies, an entire range of the same Raman frequencies can be 
excited and studied. Previously, narrowband lasers were the 
primary source available for stimulating CARS, so that the 
intensity modulation was created by the beats between two 
narrowband pulses. Broadband sources have the benefit of 
providing a much more easily tunable source of beats of a 
particular frequency, because the laser need not be retuned, 
only the pulse shaped differently. 

An ultrabroadband source of pulsed radiation can be used 
to both stimulate CARS and act as the reference. The pulses 
from the source will be divided by a frequency- selective 
element such as a dichroic beam splitter into lower and higher 
frequency pulses. The higher frequency pulse bandwidth will 
correspond to the anti-Stokes frequencies emitted by the 
sample, and will act as the reference pulse to demodulate the 
anti-Stokes signal. The lower frequency pulse will be shaped 
to stimulate the Raman frequencies of the sample in a par- 
ticular bandwidth. This scheme is illustrated in FIG. 22 . 

The strategy used here is to split the lower frequency pulse 
into two copies. Each copy will be put through separate dis- 
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persive elements, and then recombined afterwards with a time 
delay between them. Dispersive elements can be imple- 
mented by one or a combination of dispersive transparent 
materials (e.g. optical glasses or liquids), prism or grating 
pulse compressors or expanders, spatial-light-modulator 
based pulse shapers (e.g. liquid-crystal or acousto-optic Fou- 
rier plane pulse shapers), or dispersive optical fibers or 
waveguides. Dispersion causes the various frequency com- 
ponents of the pulses to be spread out in time. When the two 
pulses are overlapped, a beat frequency is produced at the 
difference between the instantaneous frequencies of the two 
pulses at a given instant. By causing the frequency difference 
of the two pulses to vary between a lower and higher differ- 
ence frequency during the time interval they are overlapped, 
the beats can stimulate the Raman frequencies in the same 
range. By varying the instantaneous frequency of the two 
pulses, and their relative delay, the range of Raman frequen- 
cies the overlapped pulses stimulate can be varied. 

To see how two dispersed pulses can produce a beat fre- 
quency spectrum that ranges from Q L <Q<Q H , consider two 
band limited pulses with different dispersion phases imposed 
on them: 

E l ((a)=E l expfKj^G))) for (o 0 -Aco/2<co<co 0 +A(o/2 (21) 

Ei(oo)=0 otherwise 
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dt t=T/2 dt t=T/2 
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d dfi 


The inequalities of Eq. 27 ensure that the beat frequency 
remains inside the excitation interval. The instantaneous fre- 
quency of two pulses that fit these criteria are shown in FIG. 
23 : 

Based on the pulses of Eqs. 21 and 22, we can determine a 
method of finding the CARS signal. Based on Eqs. 1 and 2, 
these will be: 


25 
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E 2 (m)=E 2 exp(z'(J) 2 ((j))) for (o 0 -A(o/2<co<(0 0 4A(i)/2 (22) 


E 2 (co)= 0 otherwise 

The dispersion phase (^(oo) and 4 > 2 (to) correspond to the 30 
total phase a frequency to of its respective pulse accumulates 
in its respective dispersive pulse shaper. For example, if pulse 
#1 travels through a medium with dispersion relation k(co) 
and thickness d, then the dispersion phase for that pulse 
(|q(oo)=k(oo)d . In the stationary phase approximation, the 35 
time-domain signals for these pulses are: 


The last equation can be recast as: 


E A (t) = I* 3 \t)\ 
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Ei (w) = 0 otherwise 
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dco dt 


E 2 (t) = E 2 
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exp(M> 2 (r)) for r 2 (w 0 - Aw/ 2) < 

t < r 2 (w 0 + Aw/2) 


(24) 


E t {oj) = 0 otherwise 


d(p2 d<b 2 

where — — = r 2 (w), co 2 (t) = r 2 (w), and — — = w 2 (t). 
d w dt 


We wish to further restrict the time interval of the overlap 
be confined from -T/2<t<T/2. To do this and utilize the full 
bandwidth of the signal we force pulse 1 to end at time T/2, 
and force pulse 2 to begin at time -T/2: 


The inversion in all cases including that of Eqs. 29 and 30 
can be implemented by the regularized least-squares solution 
40 already described. In general, because the anti -Stokes radia- 
tion generated by pulse #1 and pulse #2 overlap in frequency, 
a solution like that implemented for Eq. 1 5 will not suffice for 
this case. 

We note that if the lowest Raman frequency is stimu- 
45 lated exactly at the end of the overlap of the two pulses at time 
t=T/2, then there is no time to accumulate anti-Stokes signal 
from that Raman frequency. Therefore in practice one should 
choose an slightly lower than the minimum desired 
Raman frequency, perhaps decreased by 20% of 
50 This ensures that the upper frequency pulse continues long 
enough to read out the minimum desired Raman frequency. 

It is also possible to reverse the beat frequency so that the 
frequency is swept from Q L to Q H . To do this, one can pose the 
following constraint on the signals: 

55 
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0 


We would also like the frequency difference between the 
two pulses to start at and end at Q L . This best places the 65 
anti- Stokes frequencies outside of the bandwidth of the pump 
and Stokes frequencies: 


with the added constraint of Eq. 27. In addition, one can 
filter the output signal from the sample for frequencies less 
than the minimum frequency of the illumination oo 0 -Aco/2. In 
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this case, one captures the CSRS output of the sample. A 
reference pulse of this same frequency band will then 
demodulate this signal, and a regularized least-squares inver- 
sion operator for Eq. 7 can be employed to find the Raman 
spectrum. 

As an example of a pulse combination that can stimulate a 
band of Raman frequencies [2] , we design a pulse that sweeps 
the beat frequency linearly from to Q H . To do this, we will 

split a transform-limited pulse of bandwidth co 0 -Aco/ 
2<co<co 0 +Aco/2 into two components. Each component will 
have a linear chirp (quadratic phase) applied to it in the 
frequency domain. This chirp is imparted by the dispersive 
pulse shaper elements of the device. To impart a linear chirp, 
a pulse dispersing system that could consist of dispersive 
materials (e.g. optical glasses or crystals), prism compressors 
and expanders, grating compressors and expanders, disper- 
sive waveguides and optical fibers, and Fourier-plane pulse 
shapers will be needed. Producing a linear chirp will likely 
require cascaded positive and negatively dispersive elements 
to cancel the higher dispersion orders (e.g. cubic and quadric 
dispersion). The design and characterization of particular 
phase profiles 4>(oo) is well documented in the literature. 

The spectrum of the pulse combination that produces beats 
from to Q H over an interval T is given by: 



Ei(oo)=0 otherwise 

where a=(2Ao)-Q i7 -Q L )/2T, (3=(Q i7 -Q i )/2T, and t=(T/ 2) 
[ Qfy (A co - Q^) + Q L / ( A oj- £fy)] . A cosine apodization window 
has been added to improve the stability to the inverse. The 
constant a is the common chirp to both pulses, (3 is the 
difference chirp between both pulses, x is the time delay 
between the pulses, and k is the difference in field magnitude 
between the pulses. To be able to form beats at all Raman 
frequencies, Aco>Q i7 . FIG. 24 shows a graph of the frequency 
vs. time for this pulse. 

As illustrated in FIG. 2, a preferred embodiment NIVI 200 
has the following components [1]: 

Oscillator. This oscillator 201 produces an optical field that 
can excite the resonance mode of the target molecule through 
a nonlinear technique (usually stimulated Raman scattering), 
and also the photon to stimulate the photon that is measured 
(also usually through stimulated Raman scattering). The 
combination of these processes is called CARS. 

Reference Generator. The reference signal generator 203, 
which can sometimes be incorporated into the oscillator, con- 
verts part of the oscillator signal to a reference signal that can 
be used in the demodulator section. It acts as a known signal 
that demodulates the unknown signal from the sample in the 
interferometer. 

Microscope. The microscope 205 delivers the field pro- 
duced by the oscillator to the sample, and collects the field 
emitted by the sample. The excitation field is usually deliv- 
ered by a microscope objective, where the oscillator field is 
focused tightly or sparsely, depending on the scanning 
method. This focus is scanned through the tissue, and based 


on the signal received from each tissue volume an image can 
be formed. When the oscillator signal is delivered to the tissue 
containing a molecule with a compatible resonance, a non- 
linear process such as CARS can occur and produce a new 
5 sample signal (called the “anti-Stokes” for CARS processes). 
This sample signal serves as an indicator of the presence of 
the molecular resonance, and also provides additional infor- 
mation about the molecule through the temporal structure of 
the sample signal. 

to Demodulator. The demodulator 207 combines the signal 
received from the sample with the reference signal. This is 
typically achieved by constructing an interferometric cross - 
correlator. The cross -correlation of the two signals is then 
measured by a single photodetector or array of photodetec - 
15 tors. The power received by these photodetectors allows the 
cross -correlation signal to be inferred, from which the tem- 
poral response signal from the sample can be also inferred. 
With knowledge of the physics of the molecule, the presence 
of and potentially properties of the molecule being tested can 
20 be inferred from its temporal response. 

Recorder. The data recorder 209 records the data measured 
by the demodulator. This data can be digitally processed to 
produce an image that a human operator can interpret. 

Each of these modules can be implemented in a variety of 
25 different ways that can be tailored to various data acquisition 
needs. In addition, while these units are the basic units of the 
invention, often the parts can be consolidated to simplify 
implementation or reduce cost. While the basic block design 
could be construed as that of a standard interferometric 
30 microscope, nonlinear processes are occurring in the refer- 
ence generator and “Microscope” sections that allow the reso- 
nance information of the sample to be pumped. 

Each of these units will be detailed presently. 

1 . Oscillator. 

35 The oscillator produces the electromagnetic field that 
stimulates the excitation of the resonance to the probed. It 
also provides the photon that stimulates the output photon that 
is detected as evidence of CARS or CSRS. There are many 
types of oscillators and fields that can produce CARS. Each 
40 pulse produced by the oscillator should be nearly identical so 
that it can excite consistent signals in the reference generator 
and sample. If the oscillator produces too variable of a signal, 
the signals from the reference generator and sample may 
change and produce signals that can be confused with noise 
45 sources. Variability in the oscillator output is a noise source in 
itself that adds uncertainty to what the expected demodulated 
signal should be. 

The conventional way to produce CARS is to send in two 
overlapped optical pulses, one of which at frequency co 1 , the 
50 pump, and the other at oo 2 , the Stokes pulse, where 
(Oi-a) 2 =Aco, where Aoo is the resonance frequency of the 
molecules of interest. These pulses produce a beat frequency 
of Aco that manifests itself in the magnitude of the optical 
field. In linear time-independent optics, systems are sensitive 
55 only to the frequencies of the optical pulses themselves, and 
not any beats they may produce together. However, with 
sufficient intensity the intensity envelope may itself stimulate 
the molecule. By choosing two pulses that produce beats of 
this frequency, we can stimulate the molecule with two wave- 
60 lengths that the tissue is transparent to. Once the resonance is 
stimulated, another photon of frequency eo 1 (in CARS), or of 
frequency co 2 (in CSRS) stimulates the emission of a fourth 
photon, which is of frequency 2oo 1 -oo 2 for CARS, and 2oo 2 - 
coj for CSRS. 

65 Systems that can be used to produce these two frequencies 
are shown below in FIG. 3. A common configuration to pro- 
duce pulses of these two wavelengths that are overlapped in 
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time is to have a pulsed laser produce one of these wave- 
lengths, split off of some of its energy, and use this energy to 
produce a second pulse of a lower or higher frequency. In one 
configuration 300, a pump laser 301 pumps a dye laser 303, 
for example a doubled Nd: YAG pump laser at 532 nm pump- 5 
ing a tunable dye laser. In another possible configuration 302, 
a pump laser such as a Ti-sapphire oscillator pumps an optical 
parametric oscillator (a device that converts pulses to lower 
frequencies) 305. 

In yet another configuration 304, the pump laser pumps a 10 
regenerative amplifier 307, such as a Ti-sapphire regenerative 
amplifier. The regenerative amplifier then pumps an optical 
parametric amplifier (another frequency conversion device) 
309. Alternatively, as illustrated in configuration 306, the 
pulses of each wavelength are generated by two separate 15 
pump lasers, and the time overlap is maintained by a circuit 
311 that synchronizes the two sources. In another configura- 
tion 308, the pump laser pumps a continuum light generator 
313, generating broadband light which is filtered by a filter for 
the two wavelengths with group velocity dispersion correc- 20 
tion 315. 

While directly generating the two frequencies and super- 
imposing them to produce beats is the most common way to 
stimulate CARS, this method has some disadvantages for the 
method of NIVI. In CARS and CSRS, there are two types of 25 
generated signals. Resonant signals depend on the presence 
of a molecule of a particular resonance frequency to be 
present to generate the CARS signal. Another component, 
nonresonant CARS, does not require a particular frequency to 
perform conversion. Nonresonant CARS depends on the peak 30 
intensity in the signal, while the resonant component can 
build up from many beat periods and so therefore can be 
spread out in time. Because of this, it is advantageous to 
spread the CARS signals in time to reduce the nonresonant 
component. 35 

However, when the two signals are discretely generated 
and are transform-limited (are not chirped in time), the only 
way to broaden the signals in time is to reduce their band- 
width. To achieve sufficient power- spectral-density to cause 
efficient conversion, the pulses must either generated by a 40 
low-bandwidth laser, or much power will be wasted in filter- 
ing a higher bandwidth signal. Unfortunately, the range reso- 
lution in OCT, when temporal ranging is used, is inversely 
proportional to the illumination bandwidth. This requirement 
for high bandwidth conflicts with the requirement for small 45 
bandwidth for resonance specificity. It would be desirable to 
come up with an alternate oscillator configuration that would 
preserve the resonance specificity of the low bandwidth 
pulses, but actually utilize high bandwidth signals. 

Since the nonlinear excitation of the resonant molecule 50 
depends more on the beats produced than on the bandwidth 
used to produce them, it would be desirable to take a broad- 
band pulse and reshape it into a signal with the required beat 
frequency. Recent advancements have made pulsed sources 
of very large bandwidth. Some of the methods to do this are 55 
high-bandwidth Ti-sapphire oscillators, dispersion compen- 
sated mirror Ti-sapphire oscillators, double chirped-mirror 
Ti-sapphire oscillators, and continuum generation sources. 
The optical field produced by these sources can be shaped into 
a field with the beats at the required frequency. 60 

One such method that has been demonstrated in the litera- 
ture is shown in FIG. 4. A source of laser pulses from a laser 
source 401 is sent into a Fourier-plane pulse shaper 403 that 
utilizes a spatial-light-modulator (e.g. liquid crystal modula- 
tor or acousto-optic modulator). The Fourier-plane pulse 65 
shaper enables each frequency in the pulse to have its phase 
and/or amplitude altered. By applying the correct phase and 
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amplitude to each incoming frequency, the incoming signal 
can be convolved with an essentially arbitrary signal. The 
pulse shaper is set up to reshape the incoming pulse by apply- 
ing a period phase or amplitude perturbation in the Fourier 
domain with period Aoo/N, where N is a positive integer. This 
will transform a single pulse into a train of pulses that are 
separated in time by 2jt/Aoo. If only a phase perturbation is 
used, the power of the signal can be maximally preserved. The 
larger the integer N is, the longer the pulse train will be, and 
none of the bandwidth of the original pulse will be lost. 
However, most pulse shapers have a limited number of con- 
trollable frequencies, limiting the practical size of N. 

One advantage of spatial-light-modulator based pulse 
shapers is that there is typically a wide range of pulse shapes 
that can be achieved, and the spatial-light-modulator can 
often be controlled automatically by a computer. The com- 
puter can then adjust the spatial-light-modulator to achieve 
maximum signal from the sample in a feedback loop. This 
may allow automatic correction of dispersion or aberrations 
introduced by the optics of the system, and will permit more 
flexibility in probing the molecule because of the tunability of 
the pulse shapes. 

The pulse-shaper in the “Pulse-shaper type NIVI oscilla- 
tor” 400 is well described in the literature. It consists of two 
diffraction gratings, which disperse and recombine the fre- 
quencies, two lenses that focus each frequency to a point and 
recollimate each frequency, and a pulse shaper placed at the 
focal plane to separately operate on each frequency. The pulse 
shape is altered by dispersing each frequency to a separate 
angle, and then imaging each frequency to a separate point on 
the spatial -light-modulator. Alternatively, an etalon may be 
used to shape the amplitude of the pulse periodically. Unfor- 
tunately, while this would be simpler, it modifies the spectrum 
of the pulse and therefore introduces artifacts into the NIVI 
image. 

An alternative method is to take a pulse and impart a linear 
chirp to it. A linear chirp turns a pulse into one where the 
frequency rises or falls at a linear rate as a function of time. 
This rate is characterized by a constant a, which is the change 
in frequency per unit time. It is called “chirped” because of 
the noise of the equivalent sound wave. If two copies of the 
chirped pulse are created, delayed with respect to each other 
by imparting a variable time delay 407, and recombined, the 
resulting pulse will have two simultaneous frequencies that 
will rise or fall together at the same linear rate, but always be 
separated at a given instant by the same frequency. If this 
separation frequency is chosen to be Aw, then the envelope of 
the pulse will be modulated by beats of this frequency. This 
method is especially convenient because the probed reso- 
nance frequency can be adjusted easily by adjusting the time 
delay between the two chirped pulses, which is relatively easy 
and inexpensive. This will enable a NIVI instrument that can 
be rapidly and easily adjusted to scan a wide range of molecu- 
lar resonances. Systems based on tunable frequency sources 
will likely be much more difficult to dynamically change 
reliably and often. 

The ability to linearly chirp a pulse is well known in the 
literature. It can be accomplished with a pulse shaper 405 
having a combination of prisms, diffraction gratings, lenses, 
mirrors, and/or dispersive materials. Combinations may be 
required to ensure that the resulting chirp is linear and does 
not contain significant amounts of higher-order dispersion. 
Higher-order dispersion would limit the resolution to which 
the resonance could be addressed and exclude other nearby 
frequency resonances. In a typical Chirped CARS NIVI setup 
402, the chirp rate required would be fixed and the chirp rate 
should need little or no adjustment in the field. Measuring 
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devices such as Frequency Resolved Optical Gating can test 
whether a chirped pulse is linearly chirped. 

When using high-bandwidth excitation for CARS, it is 
important to filter out the entire bandwidth of excitation 
before detection so it does not interfere with detection of the 
emitted anti-Stokes light (Stokes for CSRS), because the 
nonlinear emission can not be easily distinguished from the 
much larger linearly scattered excitation light. However, this 
linearly scattered light contains the same structure that con- 
ventional OCT imaging does, and may be used to measure 
this information at the same time that a NIVI image is 
recorded. This will be convenient for superimposing OCT 
and NIVI data onto the same image, because the acquisition 
of both types of data can be designed into the same instru- 
ment. This can be implemented in practice by using a dichroic 
beamsplitter to separate the excitation and response radiation, 
and measuring the cross -correlation of the two frequency 
bands separately using a cross-correlation demodulator. 

With high-bandwidth sources where the entire bandwidth 
need not be utilized to produce the excitation field, it is 
possible to use the upper end (for anti-Stokes) or lower end 
(for Stokes) of this bandwidth as a reference field, eliminating 
the need to separately generate a reference field. However, the 
frequencies of the reference must occupy the same band as the 
received CARS/CSRS light from the sample. Some sources, 
especially continuum generation sources, will likely generate 
much more bandwidth than needed for pulse shaping and 
therefore will probably have this extra bandwidth available 
for this use. While a different process from CARS/CSRS 
typically generates this light, it will likely remain phase- 
coherent with the CARS/CSRS light and therefore should be 
useful as a reference. Phase-coherence depends on the 
mechanism of pulse/continuum generation and therefore its 
phase-coherence stability properties of a particular source 
type must be established before it is suitable for this purpose. 
A dichroic beamsplitter may be used to separate the fre- 
quency band corresponding to the response radiation from the 
oscillator energy, so that it may be utilized as a reference 
signal. 

It is also possible to simultaneously stimulate the excita- 
tion of several resonances if the sample is illuminated with the 
appropriate pump and Stokes/anti- Stokes beams. For 
example, in CARS a narrowband pump signal and a wideband 
Stokes can be used to address many resonances simulta- 
neously. This is called multiplex CARS and can be extended 
to CSRS with a broad anti-Stokes wavelength range. This 
may be used to measure the presence of several molecular 
resonances simultaneously in the sample. In addition, if sev- 
eral excitations can be produced in the same molecule simul- 
taneously, the molecule will evolve to various quantum states 
depending on the relative amplitude and timing between the 
CARS/CSRS stimulating signals for each resonance. This 
may be produced by sending in multiple pairs of Stokes/anti- 
Stokes wavelengths and pump beams in with varying time 
delays between them. By varying the time delay between 
excitations, the molecule can be made to prefer Stokes or 
anti- Stokes emissions from a particular quantum state. This 
way, the amount of anti- Stokes radiation generated from each 
quantum state could be potentially probed to identity the 
molecule. 

In general, the spatial -light-modulator system of FIG. 4 
could be used to produce more general pulse shapes than a 
series of beats at a single resonance frequency. By using a 
more complicated pulse shape, several bonds present in a 
molecule can be coherently excited, and energy transferred 
between them in a coherent fashion. Because each molecule 
has some difference between the bonds present and their 


18 

relative orientation (and therefore the coupling in the quan- 
tum wave functions between them), pulses can be shaped that 
will selectively transfer energy between the states for a par- 
ticular molecule, and not be selective for other similar mol- 
5 ecules. In this way, the emission of stimulated Raman scat- 
tering or another coherent scattering process can be made 
more specific than just every molecule possessing a bond of a 
particular energy. With an automatically controlled pulse 
shaper, such as those based on spatial light modulators, feed- 
10 back can be employed where the computer can test various 
pulse shapes, measure the resulting emitted light temporal 
signal using the demodulator, and progressive reshape the 
pulse to optimize the signal from the target molecule and 
exclude other molecules. Once a usefi.il temporal field shape 
1 5 for stimulating a molecule has been found, it can be stored in 
a database and later used for identifying that molecule in the 
future. 

2. Reference Generator. 

20 The reference generator takes a portion of the signal pro- 
duced by the oscillator and converts it to a reference signal. 
This reference signal is later used to demodulate the sample 
signal. The reference generator is a nonlinear process that 
converts light in the illumination bandwidth to light in the 
25 sample’s emission bandwidth, so that interference can occur 
between them. This nonlinear process may or may not be 
CARS or CSRS. 

A common implementation of the Reference Generator 
would be to focus the oscillator excitation into a sample of the 
30 same molecular species that one wishes to image. The refer- 
ence signal should then be very similar to the same molecular 
species contained in the sample. This is because they are the 
same molecule, illuminated by nearly identical pulses, con- 
verting them to the output signal using resonant CARS or 
35 CSRS. The magnitude of the cross-correlation between these 
two signals should be great because of their similarity. In 
addition, if there is variability of oscillator output, having the 
reference generator and sample contain the same substance 
will respond in similar ways, so that the cross-correlation 
40 signal can remain high despite fluctuations in the oscillator. 
The benefit of using the same molecule in the reference gen- 
erator is that it is the molecule’s signal that acts as its own 
“fingerprint” with which the cross-correlation can use to rec- 
ognize the molecule in the sample. If more selective excita- 
45 tion processes than CARS are used, then using the same 
molecule in both reference and sample will help ensure that a 
reference signal can be generated for a given excitation field. 

A nonresonant nonlinearity can be used as the reference 
generator as long as the peak power of the oscillator signal 
50 can excite a sufficient quantity of reference signal . High peak 
power can be maintained by not chirping the oscillator signal 
that is sent to the reference generator, while sending a rela- 
tively low peak power signal to the sample. Nonresonant 
CARS or CSRS can be used to generate an anti-Stokes or 
55 Stokes signal, respectively, from a medium that does not 
necessarily have a resonance at the frequency of the target 
molecule. The benefit of this is that the medium may not have 
to be changed each time a different molecular species is 
scanned for, because otherw ise a medium with a resonance at 
60 that wavelength would need to be chosen. Also, this species 
can act as a standard signal source against which the return 
signals from many samples can be compared. The nonreso- 
nant CARS can be implemented by focusing the excitation 
radiation into a sample of liquid that produces a CARS/CSRS 
65 signal in the same frequency band as that generated from the 
sample. For example, benzene will generate a CARS anti- 
Stokes signal in the 3000-3100 cm -1 frequency band. 
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Continuum generation is another type of nonresonant non- 
linear process that can be used in the reference generator. A 
sufficiently high peak power pulse is focused into a medium, 
where it excites a broad bandwidth of frequencies to be pro- 
duced. If the produced frequencies overlap the emission fre- 5 
quency band produced in the sample, this portion of the 
continuum can act as a reference signal. The generated con- 
tinuum must be created by a mechanism that is sufficiently 
stable to not be overly sensitive to fluctuations in oscillator 
intensity. An unstable reference signal will result in noise in 10 
the cross-correlation signal. The benefit of continuum gen- 
eration is that is likely to create a broad bandwidth that will 
have signal in the emission bandwidth of the sample, so that 
the continuum need only be filtered for the needed frequency 
band. Also, if the oscillator employs continuum generation, it 15 
may already generate light within the emission bandwidth 
that can be used as a reference generator, eliminating a sepa- 
rate nonlinear process in the reference generator step. Some 
examples of materials used for continuum generation mate- 
rials are optical glass, fused silica, calcium fluoride, sapphire, 20 
ethylene glycol, water, high numerical aperture optical fibers, 
photonic crystal optical fibers, microstructured optical fibers, 
dispersion-shifted optical fibers, and gas cells (e.g. cells filled 
with helium, argon, or nitrogen). 

Other candidate processes for nonresonant nonlinear ref- 25 
erence generation include second and higher harmonic gen- 
eration, stimulated Raman scattering, sum and difference 
frequency generation, optical parametric amplification, four- 
wave mixing, and self phase modulation. 

FIG. 5 shows an implementation of a reference signal 30 
generator 500. The concentration optics 503 are typically 
implemented as some combination of lenses and mirrors. The 
concentration optics may also require some combination of 
frequency dispersive elements such as prisms, diffraction 
gratings, pulse shapers, and dispersive materials to prepare 35 
the temporal shape of the signal for nonlinear generation. 
Concentration in space and time may be necessary because 
the nonlinear processes are power sensitive, and depends on 
the strength of the nonlinear process. When the light 501 
enters the nonlinear medium 505, it undergoes conversion to 40 
a frequency band coinciding with the frequency band of the 
response signal from the sample. This nonlinear medium may 
be one of the media mentioned above, either a sample of a 
target molecule, a solvent, or a continuum generation 
medium. After exiting the nonlinear medium, the reference 45 
signal 509 is collected by reference collection relay optics 
507, where it is sent to the demodulator where it is combined 
with the sample signal. The collection relay optics are usually 
implemented as some combination of lenses and mirrors that 
collimate the reference field radiation. This reference field 50 
should be characterized to find its temporal structure by 
instruments such as Frequency Resolved Optical Gating, 
cross-correlation with another known signal, nonlinear sono- 
grams, or nonlinear autocorrelations/cross -correlations. 

3. Microscope. 55 

The microscope delivers the excitation radiation from the 
oscillator to the sample and collects the resulting sample 
emission. Inside the sample, a coherent nonlinear process 
such as CARS or CSRS takes place that, in the presence of a 
molecule of interest, will emit the sample field in response to 60 
the excitation field. The sample field is collected by the micro- 
scope and then propagated to the demodulator, where the 
sample field can be estimated from the measured cross-cor- 
relation between the sample field and reference field. 

Microscope systems can be differentiated by various 65 
implementation choices. They can either illuminate one (se- 
rial scanning) or many points (full field imaging) at a time on 
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the tissue. If the pump and Stokes (anti-Stokes) beams of the 
excitation field are separated in frequency, they can be sent in 
either separate (non-col linear) or identical (collinear) angles 
into the sample. The temporal delay of the response radiation 
relative to the reference may or may not be used to range 
molecular constituents in the tissue. The response radiation 
can be collected in the forward scattering (forward CARS/ 
CSRS) or backward scattering (epi-CARS/CSRS) directions. 

The microscope measures spatially resolved molecular 
density by illuminating various points on the tissue with the 
oscillator field, collecting the emitted sample field, and 
recombining with the reference field in the demodulator. One 
point in the tissue may be illuminated at a given time, result- 
ing in serial or raster scanning of the molecular density 
through the tissue. Alternatively, a line or a complete plane of 
points may be illuminated, so that data may be acquired from 
many points in parallel. Illuminating and measuring the radia- 
tion from an entire plane of points is called full field imaging. 
At the time of this writing, full field imaging is seldom used 
because it requires an array intensity detector such as a 
charge-coupled-device (CCD) to simultaneously measure the 
demodulated signals of all of the illuminated points. Unfor- 
tunately, as of this writing CCD arrays produce thermal dark 
noise at each pixel, and also have a relatively limited dynamic 
range of measurable intensity values. Demodulated interfer- 
ence signals often require very high dynamic range detection. 
It is conceivable that future CCD or other types of focal plane 
arrays (e.g. CMOS arrays) may overcome these limitations. 
Full field imaging also requires that the tissue be illuminated 
by larger amounts of power because measurable signal must 
be produced for an entire area rather than just one point. Since 
this is more likely to result in tissue damage, full field imaging 
will probably be used when speed of acquisition is para- 
mount. 

FIG. 6 shows three examples of full field CARS configu- 
rations. The “beam delivery optics” are usually implemented 
as some combination of mirrors and lenses that deliver beams 
that illuminate a wide area or line on the sample. For all of 
these microscope configurations, beam delivery and collec- 
tion optics will typically utilize a microscope objective. The 
“response field collection optics” 601 are similarly imple- 
mented as a combination of lenses and mirrors that relay the 
response field to the demodulator so that it may be recom- 
bined with the reference and detected. The noncollinear full 
field CARS 602 delivers the pump 603 and the Stokes 605 
beams (assuming they are separate) at separate angles, so that 
the response is separated in angle from the illumination by an 
angle given by phase-matching conditions. The collinear 
geometry 607 sends the pump and Stokes radiation in the 
same directions (or in a single field if they can not be sepa- 
rated) and collects the radiation in the same direction, which 
can be discriminated with a dichroic beamsplitter. The epi- 
CARS geometry 611 collects the backscattered radiation, 
usually through the same objective optics that the sample is 
illuminated through. The epi-CARS can be discriminated 
from the illumination with a dichroic beamsplitter or an inter- 
ference filter. 

If collinear CARS is used, where the pump and Stokes 
(anti-Stokes) beams from the excitation beams overlap, then 
the response field 609 will overlap the excitation beams. Then 
the response field frequency band should not completely 
overlap the excitation frequency band, so that a spectral filter 
may distinguish between the excitation and response fields. If 
non-collinear CARS is used, then the response beam can be 
sufficiently angularly separated (as determined by the phase- 
matching criterion) from the excitation radiation to be filtered 
by a spatial filter. However, in the case of non-phase-matched 
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CARS, such as epi-CARS, the interaction CARS volume 
must be small enough so that the response is radiated effec- 
tively isotropically, so that spatial filtering is unavailable and 
spectral filtering should be used. 

The other, more commonly used alternative is serial point 
scanning. Serial scanning tightly focuses the oscillation sig- 
nal into the tissue to create a very small volume where peak 
power is maximized. The focusing is usually achieved using 
a microscope objective. This focus is then scanned through a 
3-D set of points in the tissue, and the sample signal gathered 
from each point is demodulated to produce a 3-D NIVI image. 
Since nonlinear processes are power sensitive, efficient 
CARS/CSRS occurs only at the focus. If the sample is small 
enough, the focus may be scanned through the sample by 
translating the sample in all three dimensions. However, it is 
not feasible to move large samples such as human subjects 
this way. In this case, the beam focus can be moved in the 
transverse direction by steering the beam, perhaps using gal- 
vanometer rotated mirrors, acousto-optic modulators, or 
translating the lens assembly. The depth may be scanned by 
mechanically adjusting the distance between the lens and the 
tissue, perhaps using a lead-screw translator and/or a piezo- 
electric transducer. Since the signal can be excited at only one 
point at a time, one can be sure that the resulting measured 
sample signal is due to the interference of emissions of radia- 
tion produced in that volume only. This can be a benefit in 
NIVI when high phase resolution is required because one is 
assured that any measured phase shifts are not due to inter- 
ference between molecules at disparate spatial locations. This 
higher phase resolution may be used to better differentiate 
between similar molecular species. Serial point scanning 
typically utilizes a single photodetector or a small number of 
photodetectors at the demodulator, which has the benefit that 
the dark current of a single photodetector is usually less than 
that of an entire CCD array, and a single photodetector can 
also typically handle a higher dynamic range of measure- 
ments. 

FIG. 7 shows the geometry of translated serial-port scan- 
ning configurations. In all of these configurations, the sample 
is on a translator that moves the sample through the focus to 
form an image of the molecular density at various points. The 
translator could be a three-axis linear screw drive translator, 
or piezoelectric translator, or a combination of these. The 
beam delivery optics 701 focuses the pump 703 and Stokes 
(anti-Stokes) 705 beams onto the tissue at a point of interest, 
and the response field collection optics 707 gathers the gen- 
erated anti-Stokes (Stokes) field 709 from the tissue. In the 
non-collinear geometry 700, the pump and Stokes illuminate 
the point of interest at different angles, so that the anti-Stokes 
will emerge at a third angle given by phase-matching, so that 
the anti- Stokes is spatially separated from the illumination. In 
the collinear geometry 702, the illumination and response 
fields emerge overlapped, so that they must be discriminated 
by frequency (e.g. using a dichroic beamsplitter). Finally, in 
the epi-CARS geometry 704, the backscattering response 
fields are collected, often by the same optics through which 
the illumination is projected, and can be separated with a 
dichroic beamsplitter. 

Alternatively, the focus can be moved and the sample can 
be left stationary. This can be accomplished by a combination 
of tilting the illumination beams before it enters the beam 
delivery optics and/or translating the beam delivery optics 
around (FIG. 8). Translating the collection optics and/or tilt- 
ing the exiting response beam with a beam steering optics 801 
can capture the exiting response field. The beam steering 
optics changes the direction of the incoming beam 803. This 
can be implemented, for example, by a galvanometer- 
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scanned rotating mirror. By changing the direction of the 
beam before it enters the beam delivery optics 809, the posi- 
tion of the focus in the sample can be changed. The translation 
for the beam delivery and collection optics 807 can be imple- 
5 mented with a piezoelectric and/or a linear screw-drive trans- 
lation stage. Translating the beam delivery and/or collection 
optics moves the focus with the optics through the sample. 
These two mechanisms can be combined to enable the three- 
dimensional translation of the focus through the sample. This 
to configuration is especially convenient for the epi-CARS con- 
figuration, because the same beam steering optics and beam 
delivery optics 809 can be used to collect the response field 
805. For these configurations, a compensating delay may 
need to be incorporated into the demodulator because steer- 
15 ing the beam and/or translating the delivery and collection 
optics change the travel time for the illumination and 
response signal through the microscope section. 

All of the previous scanning modes, full field imaging, 
translated serial-point scanning, and beam-steering serial- 
20 port scanning, use the spatial location of the illumination 
beam to differentiate between response signals gathered from 
various points in the sample. This confinement method is the 
same as that used by multiphoton microscopy or CARS 
microscopy. Other technologies such as Optical Coherence 
25 Tomography and Optical Coherence Microscopy use tempo- 
ral ranging in addition to spatial confinement to further isolate 
the contributions of signal from various points in the sample. 
Because of the heterodyne nature of NIVI, this scanning 
mode is also available. It may be attractive for in vivo imaging 
30 because it will enable scanning tissue without translating the 
microscope objective or sample, and may scan faster because 
there are mechanisms for scanning the temporal delay much 
faster than translating objective optics. The phase measure- 
ment capability of NIVI is useful for both temporal ranging 
35 and molecular species identification this way. 

FIG. 9 shows two configurations of NIVI 900 and 902 
utilizing temporal ranging. Temporal ranging is achieved by 
measuring the interference of the reference and response 
signals for various relative temporal delays. The temporal 
40 gating configuration superficially resembles the other epi- 
CARS configurations. However, unlike previous scanning 
modes, the depth-of-field of the focusing of the illumination 
from the beam delivery optics in the tissue is set to be long, 
because the temporal gating will be used to discriminate 
45 between molecular constituents within the depth-of-field. 
This temporal ranging utilizes epi-CARS/epi-CSRS because 
the backscattered response signal is collected from the 
sample by the response field connector so that the signal delay 
into the tissue can be timed. Because CARS/CSRS is not 
50 phase-matched for the backwards direction, the generation of 
a backscattered signal will only occur efficiently for small 
particles or edges of particles with a compatible resonance. 
The chief difference between this configuration and the other 
epi-CARS configurations is that the interference signal is 
55 scanned as a function of relative delay, and an interference 
signal maximum indicates the presence of a molecular spe- 
cies at a particular depth in the medium (corresponding to that 
time delay). In the lull field configuration 900, an entire plane 
of points is interfered with the reference signal to produce 
60 simultaneous measurements of the cross-correlation signal 
for the entire plane 909. By scanning the delay, the molecular 
density of various planes can be measured. For the beam- 
steered temporal gating setup 902, the beam is scanned 
through the sample by tilting the beam through the beam 
65 steering optics 903. Beam steering provides lateral displace- 
ment of the beam, and temporal gating provides depth infor- 
mation so that a three-dimensional volume is scanned. The 
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sample can be translated laterally also to scan the beam. The 
implementation differences between the temporal gating con- 
figurations are in the choice of depth-of-field of the illumina- 
tion and collection optics, and in the demodulator. The 
demodulator must be designed to scan a sufficient temporal 
interval to capture the interference signal between the refer- 
ence and response signals. This temporal interval is typically 
between 500 microns to 10 mm. 

4. Demodulator. 

The purpose of the demodulator section is to decode the 
response signal from the sample so that it may be measured by 
relatively slow electronic equipment (slow compared to the 
oscillations of the electric field of the response signal). The 
magnitude of this demodulation signal will be related to the 
density of a molecular species of interest in the tissue. With 
knowledge of the molecular density at each point, a molecular 
density map, or NIVI image, can be presented to the user. This 
demodulation is implemented as the cross -correlation of the 
response signal relative to a generated reference signal. Out- 
lined below are various optical configurations that produce 
this cross-correlation signal. 

There are various design choices that are made when 
choosing a demodulator. First, one needs to know whether or 
not full field microscopy is used. Also, one must decide 
whether the cross-correlation will be measured one sample at 
a time, or many samples in a single instant. 

FIG. 10 shows an example of a cross -correlator that can be 
used to demodulate the full field CARS signal, measured for 
example as in FIG. 6. The response field is collected and 
relayed by the “response field collection optics” of FIG. 6 to 
the Response Field 1001 of FIG. 10. The reference field 1003 
is produced by the Reference Generator and relayed to the 
cross -correlator. The reference is delayed relative to the 
response field, and they are mutually overlapped using a 
beamsplitter 1005. The “response field spectral filter” 1007 
filters the recombined field for only the frequency band that 
contains the response field bandwidth, which removes any 
remaining oscillator signal. The filter can be implemented by 
a combination of interference filters, color glass filters, dich- 
roic beamsplitters, or other frequency selective elements. The 
combined fields are imaged onto a focal plane array 1011, 
such as a CCD, so that the CCD is the conjugate image plane 
of the sample plane. The imaging is achieved with the “relay 
imaging optics” 1009 which are some combination of lenses 
and mirrors. The intensity detected on the CCD corresponds 
to a single cross -correlation signal measurement collected 
from various points in the sample. The “variable attenuator” 
1013 is adjusted to maximize the use of the dynamic range of 
intensity measurements of the CCD. As the focus in the 
sample changes (e.g. by translating the sample), and/or the 
adjustable delay 1015 is changed, the entire cross -correlation 
signal can be measured, forming a three-dimensional data set. 
The typical scan range length for high numerical aperture full 
field imaging will be up to 100 microns. The delay can be 
implemented as a mirror translated by a linear-screw-drive 
translation stage, or by piezoelectric actuation. This same 
cross -correlation can also be used for temporal ranging full 
field NIVI, as shown in FIG. 9, with the only difference being 
that the delay mechanism must be designed to scan a suffi- 
ciently long range of interest in the sample, typically from 
500 microns to 10 mm. 

The measurement of a cross-correlation signal for a serial- 
point scanning microscope is basically the same system as for 
full field imaging, except that a single photodetector can be 
used rather than an array photodetector. The adjustable delay, 
response field spectral filter, and relay imaging optics can all 
be implemented in similar ways to the full field case. The 
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relay imaging optics need only focus the combined response/ 
reference signal onto the photodetector. In some cases, when 
only the magnitude of the cross-correlation signal at its peak 
needs to be measured, it may be desirable to dither the adjust - 
5 able delay with a piezoelectric transducer a fraction of a 
wavelength, so that the magnitude of the cross -correlation 
signal peak can be demodulated with a multiplying mixer and 
low-pass filter (an electronic heterodyne demodulator). This 
configuration will be covered in more detail in part five. To 
to measure the entire cross -correlation, the adjustable delay will 
be scanned over various time delays and the photodetector 
intensity signal measured. 

In FIGS. 10 and 11, a variable attenuator 1013 is used to 
adjust the intensity of the signal that reaches the photodetec - 
15 tor 1101 so that its dynamic range is not exceeded. Variable 
attenuators can be implemented with liquid crystal shutters, 
neutral density filter wheels, or rotating polarizers. The spa- 
tial filter 1103 is used on the response signal to filter out 
spatial inhomogeneities in the response signal that could 
20 reduce the depth of modulation at the photodetector. A spatial 
filter would typically consist of a telescope of two converging 
lenses, with a pinhole in the focal plane between the lenses to 
filter out power around the main diffraction focus. 

When the oscillator produces pulses of a low repetition 
25 rate, so that for example high peak power can be employed, 
one may want to measure multiple samples of the cross- 
correlation signal simultaneously. This can be achieved by the 
configurations in FIG. 12. These configurations have the 
advantage that, with a sufficient amount of pulse power and 
30 number of photodetector array samples, the cross -correlation 
signal can be measured using a single pulse. Low repetition 
rate illumination can keep the peak power high while average 
power remains low. 

The “linear photodetector array cross-correlator” 1201 
35 expands the response and reference signals, and interferes 
them with an angle between the two beams. The beam expan- 
sion is achieved by, for example, a pair of converging achro- 
matic lenses 1203 arranged in a telescope configuration. The 
points on the wavefront where the two signals combine will 
40 have various time delays between them. A cylindrical lens 
1205 then focuses the beams into a line image on a linear 
photodiode array 1207. This concentrates the signal to adjust 
for the narrow height of the detector pixels. Each intensity 
sample on the linear photodetector array indicates a sample of 
45 cross -correlation of the two fields with various relative time 
delays, with a constant intensity signal added. The recorder 
would then read the linear CCD signal so that a complete 
NIVI image can be assembled. 

The primary disadvantage of this system is that the modu- 
50 lation of the intensity signal on the linear photodetector array 
is rather low, and the dynamic range of typical CCDs is 
likewise low, so that the signal cannot be measured to a high 
signal-to -noise ratio. To combat this, one can utilize the Fou- 
rier-transform cross-correlator 1209. Rather than directly 
55 measuring samples of the cross-correlation, it interferes the 
two signals and measures their spectral decomposition. The 
beamsplitter 1005 combines the response and reference sig- 
nals, which are filtered for the response bandwidth by utiliz- 
ing a response field spectral filter 1007. This signal is then 
60 filtered by a frequency dispersive element 1211 such as a 
diffraction grating that scatters each frequency to a different 
angle. The power of each frequency is then focused to a 
different pixel on the linear photodetector array 1215 by using 
a focusing element 1213 (typically a combination of lenses or 
65 mirrors). The samples of the intensity measured on the pho- 
todetector indicate the real part of the Fourier transform of the 
cross -correlation signal. The recorder may compute an 
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inverse discrete Fourier transform of the linear photodetector 
array intensity samples to recover the cross-correlation. The 
Fourier transform already performs the necessary Hilbert 
transform to infer the imaginary part of the cross -correlation 
signal. However, because it is an intensity measurement, the 
cross -correlation signal, a time-reversed version, and the 
autocorrelation of the response signal are superimposed in the 
intensity signal. By choosing the adjustable delay longer than 
the length of the cross-correlation temporal signal, the recon- 
struction of these signal components will not overlap in the 
time domain, and so the recorder device may distinguish the 
cross -correlation from its mirror image and the autocorrela- 
tion. 

The adjustable delay may also be dithered a small amount 
to introduce a phase shift into each measured frequency com- 
ponent at the linear photodetector array. By utilizing at least 
three linearly independent phase shifts, the phase of each 
Fourier component can be known and the time domain cross- 
correlation computed using an inverse discrete Fourier trans- 
form. However, then the ability to measure the cross -correla- 
tion of a single pulse will be lost. Another possibility is to use 
a linear photodetector array with several rows (at least three) 
rather than just one row of pixels. By tilting one of the wave- 
fronts slightly vertically with respect to the other, a small 
phase shift can be introduced in the measured cross -correla- 
tion signal between rows on the linear photodetector. With a 
sufficiently large phase shift, the complex amplitude of each 
Fourier component can be inferred from the intensity samples 
from each column on the linear photodetector by utilizing the 
discrete Fourier transform. Since all of the time delayed 
cross -correlation samples are measured simultaneously, this 
is a single-pulse measurement. 

To minimize the effect of dark charge built up in the linear 
CCD detector, it is best to discard the charge as soon as 
possible before the pulse arrives, and read the CCD as soon as 
possible after the pulse arrives. The speed of readout should 
be as fast as possible given limitations in accuracy due to the 
readout noise. 

5. Recorder. 

The recorder accumulates the samples of the cross-corre- 
lation signal gathered from various points in the sample, 
processes this data, and presents these as a human-interpret- 
able image. It is usually implemented as a data acquisition 
digital computer with some method of automatic control of 
the adjustable mechanisms (such as delay lines, piezoelectric 
actuators, and galvanometer mirrors), analog-to -digital con- 
version, and some form of image output device such as a 
screen or printer. It may also have control of the oscillator 
itself, to automatically tune the wavelengths or bandwidths of 
output or control the rate or timing of pulse output. It may also 
control the delay lines or spatial light modulators in the pulse 
shaping mechanisms of FIG. 4. 

Typically, the recorder will scan the illumination through 
the sample and/or the adjustable time delay and measure the 
cross-correlation signal. Unless the Fourier-Transfoim cross- 
correlator is used, the intensity samples represent the cross- 
correlation signal with a constant level added that could be 
discarded. The magnitude of the cross -correlation signal indi- 
cates the presence of the target molecular species. For CARS/ 
CSRS the magnitude of the cross -correlation signal is related 
to the second power of the molecular density of that species. 
By demodulating the magnitude of the cross -correlation sig- 
nal, the molecular density can be estimated for the points on 
the cross -correlation signal for the areas from which the 
response signal was collected to form that cross -correlation 
signal. Both spatial confinement and temporal ranging can be 
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used to differentiate between the signals due to molecular 
densities at different locations in the sample. 

FIG. 13 documents the relationship between the elements 
of the recorder 1300. For serial -scanning or temporal-scan- 
5 ning configurations, the instructions from the operator are 
entered via a human interface device 1304 into the digital 
computer 1301. The digital computer controls the scanned 
beam position in the sample using the galvanometer mirror 
angles of the galvanometer scanned mirrors 1303, the oscil- 
lator 1302, the adjustable delay line 1309, the position of the 
microscope objective position 1305, and/or the position of the 
sample translator 1307 to scan the illumination through the 
sample. In general only a subset of these need be controlled to 
15 scan the three-dimensional volume of the tissue. If full field 
imaging is used, usually only the depth and/or the delay line 
need be scanned. The cross-correlation signal as measured by 
the intensity is converted to a voltage by the photodetector 
1311, which is in turn converted to a digital sample by the 
20 analog-to-digital converter 1313. The digital computer varies 
the delay line and/or reads various pixels from the photode- 
tector (if it is an array detector) to determine the cross-corre- 
lation signal. If a Fourier-transform cross-correlator is used, 
the computer will need to compute the inverse discrete Fou- 
rier transform of the signal to find its cross-correlation from 
the intensity samples. The computer then associates the mag- 
nitude of the cross-correlation signal with the molecular den- 
sity, and assembles these magnitudes into a density map of the 
sample. This density map is stored in the storage device 1312 
30 presented on the visual display 1315 and/or made into a 
physical representation with a printer or stereolithography 
device 1317. 

To aid in the measurement of the magnitude of the cross - 
35 correlation signal, the configuration 1400 of FIG. 14 is sug- 
gested. 

Because the digital computer 1301 can record only rela- 
tively slow signals, a dither oscillator 1401 may be used to add 
a small high-frequency dither signal into the adjustable delay 
40 line 1309 that produces a periodic perturbation of the delay in 
the signal of a magnitude usually of less than one wavelength. 
This same signal is multiplied by the received photodetector 
1403 interference signal to demodulate it, and is low pass 
filtered through the low pass filters 1405 to remove harmonics 
45 of the dither frequency. Both the dither and its quadrature 
signal are demodulated, because these correspond to the real 
and imaginary parts of the complex cross-correlation ampli- 
tude. The analog-to-digital converter 1313 will then be 
directly measuring a quantity proportional to the complex 
50 magnitude of the cross -correlation for the delay line position. 
The computer can utilize the digitized real and imaginary 
cross -correlation components to display the amplitude and 
phase of the cross -correlation signal. The amplitude will cor- 
respond to the magnitude of the reflection, and the phase will 
55 correspond to the Doppler shift of the reflection. If a very fast 
dither signal is desired (above 10 kHz), an electro-optic 
modulator or acoustical-optic modulator can be placed in the 
adjustable delay line system to rapidly modulate the delay a 
small amount. 

60 The magnitude of the cross-correlation signal depends 
both on the density of molecules available to produce the 
response signal, and on the temporal structure of the signals 
radiated from the molecules. The temporal signal produced 
by the molecules can be inferred from the cross -correlation 
65 signal. If f(t) is the temporal response signal produced by a 
molecule, and g(t) is the known reference signal, the mea- 
sured cross-correlation T(t) is given by: 
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£°f (r)g(t 


- r)dt 
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If F(oo), G(oo), and r(co) are the Fourier transforms of f(t), 
g(t), and T(x) respectively, then f (co)=F(cl))G(-cl)). The func- 
tion T(oo) can be computed from the Fourier transform of the 
cross-correlation and G(co) can be computed from the mea- io 
sured reference signal. The Fourier transform F(co) of tempo- 
ral signal f(t) can then be estimated by F(co)=r(co)G(-co)*/( 
G(-oo) 2 +N(oo) 2 ) (change to equation), where N(co) is an esti- 
mated power- spectral-density of the noise. An inverse Fourier 
transform of F(oo) yields f(t). 15 

The recorder can perform this computation and thus 
recover f(t) for various molecular species. An unknown mol- 
ecule may be identified by comparing its temporal signal to 
known molecular signals. In addition, by looking at the fre- 
quencies co of phase arg F(oo) where the phase changed rap- 20 
idly, the resonance frequencies of the molecules of interest 
can be determined to high precision. A library of the temporal 
responses and resonance frequencies of various molecules for 
CARS excitation may be built up and used to identify 
unknown molecules in vivo. We also note that the same cross- 25 
correlation measurement configurations can be used to mea- 
sure the OCT backscattered signal due to the scattering of the 
excitation radiation. To do this, the excitation and response 
bands should be separated with a dichroic beamsplitter, and 
the cross-correlation setups implemented separately to elimi- 30 
nate photon noise from the excitation band being measured in 
the response band. Different cross-correlation setups may be 
used for each configuration, e.g. a single-pulse measurement 
system for the conventional excitation and a time delay cross- 
correlator as in FIG. 11. 35 

EXAMPLES 

Theory of Cars Interferometry and Inference 

We consider a focused optical pulse, with an electric field 40 
at the focus given in the frequency domain by E z (co). At the 
focus there is a medium with a complex Raman spectrum 
given by y u (3) (Q) where Q is the Raman frequency. The anti- 
Stokes/Stokes electric field produced by the CARS/CSRS 
(Coherent Stokes Raman Scattering) process at the focus is 
given by: 


/* 3) (ft) = x (3) (^)J°° Ei (to + ft) Ei(oj)doj 

(1) 

50 

E a (to) = & 3 \£L) Ei(a) - £L)d£l (CARS) 

(2) 

E s (o>) = P®\£l)*Ei((o + ft)cfft (CSRS) 

(3) 

55 


These equations assume that classical electromagnetic 
field approximations hold, a perturbative interaction with the 
sample, and that there are no levels directly resonant with any 60 
of the individual frequencies in the pulse. These equations 
give the time evolution of a CARS/CSRS process involving 
many possible simultaneous Raman-active vibrations. In 
CARS, the molecule is excited by two frequencies and co 2 

that are separated by the Raman vibrational frequency 65 
<jOi-<jo 2 =Q. The upper frequency co 1 is called the pump, and 
the lower frequency oo 2 the Stokes. In general, the signal 
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consists of not just two discrete frequencies but a continuous 
spectrum. In this case, Eq. 1 represents the excitation of the 
molecule through SRS. A vibration at frequency Q is excited 
by every pair of optical frequencies with a frequency differ- 
ence of Q. The excitation at any given Raman frequency is 
given by the polarization P (3) (Q). This excitation is converted 
to anti-Stokes radiation E^(co) by a second SRS process 
whereby the pump signal at frequency co 1 is partially con- 
verted to the anti-Stokes signal at frequency co^=co 1 +Q=2co 1 - 
oo 2 . For broadband signals, the anti-Stokes signal must be 
summed over all pump frequencies, and so Eq. 2 results. The 
CSRS process is likewise modeled by Eq. 1 and Eq. 3. In 
general, both CARS and CSRS are generated simultaneously 
for any given pulse. 

The information about the sample is contained in its 
Raman spectrum /( 3 )(Q). By sending in a pulse with a known 
electric field E z (co), and measuring the returned anti-Stokes 
signal E^(oo), we would like to infer % (3) (Q). Unfortunately, at 
optical frequencies the electric field E^(oo) is not directly 
measurable because it oscillates on a time scale too fast to be 
electronically demodulated, necessitating interferometric 
demodulation. 

To demodulate the anti -Stokes field, we generate a refer- 
ence field R(tt)) from the oscillator field that contains frequen- 
cies same as the anti-Stokes radiation. The interferometric 
demodulator measures the intensity of the cross -correlation 
of the anti-Stokes and reference fields, which is given by the 
frequency spectrum I(oo)=R(oo) *E^(co). A good method of 
interferometric demodulation for the techniques disclosed is 
spectral interferometry [1]. Spectral interferometry may be 
especially attractive because it allows the signal to be sampled 
in one shot to minimize transient effects, and can achieve 
higher signal-to -noise ratio than conventional temporal inter- 
ferometry. Putting the steps of generating the anti-Stokes 
radiation and interferometry together, the cross -correlation 
spectrum is given by: 


7 (<u) = R / 3 Hn) £,(&> - si) 

du'dtl (CARS) 


l(a>) = R M* jH. x {3 WEi((o + ft) 

E i (oj')daj' dfl (CSRS) 


ijT Ei(oj’ + £1) Eiioj’j* 

if Ei(a)'+nf 
Jo 


(4) 


(5) 


An important consequence of Eqs. 4 and 5 is that the 
cross -correlation spectrum I(co) is a linear function of the 
Raman spectrum % (3) (Q). Because of this, linear estimation 
can be used to estimate % (3) (Q) from I(oo). To better see this, 
we reform Eqs. 4 and 5 as linear operators A and S with 
kernels A(co,Q) and S(co,Q) respectively: 


l(a>) = A X = * ,3| (fl )A(<0, Sl)dn 

where 

A (to, ft) = R (ajyEiio/ - ft)^°° Ei(co' + ft) Ei(oj')* dco' (CARS) 


(6) 
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-continued 

M =sx' = £" at 13 Ws (a, n)dn 


i 

where 


(7) 


S {co, ft) = R (coyEiico' 4 - n)J^ Ei(co' + £l)*Ei(a/)d(o' ( CSRS ) 


With the CARS interferometry process now expressed as a 
linear operator A mapping % (3) (Q) to I(co), an inverse operator 
can be found to estimate %^ 3 \Q) from I(co). Many inverse 
operators are possible and should be chosen on the basis of 
stability, ability to model the physical system and noise, and 
ability to incorporate confidence measures of the data. We 
present an inverse operator based on the weighted Tikhonov 
regularized least-squares solution. The least-squares solution 
by itself attempts to find a solution for y} 3 \Q) that minimizes 
the Euclidean distance ||I-A%|| for the linear system I=A/. 
Unfortunately, this in practice tends to produce poor results 
due to poor conditioning and the inability to incorporate 
confidence information about the data. A solution is to aug- 
ment the linear system with a weighting operator 

WI=W{(m)I{(m) ( 8 ) 

where W(oo) is a weighting function that is of unit value 
inside the anti-Stokes bandwidth and zero outside. The 
weighted linear system becomes WI=WA%. The Tikhonov 
regularized solution to this system is %=(A i7 W i AV r A+eI)“ 
iA^W vVI where the operator I is the identity operator (for 
CSRS it is x=(S i7 W // WS+6l)- 1 S i7 W i7 WI). The constant e>0 
is the set to account for the magnitude of noise in the mea- 
surement, e.g. from thermal or photon noise. In practice, this 
operator can be computed numerically using iterative linear 
solution methods such as the preconditioned conjugate gra- 
dient method. 

To perform this inverse in practice, we used the precondi- 
tioned conjugate gradient method. This method is a method of 
computing the inverse of a linear system for a particular 
solution without explicitly computing the inverse of the 
matrix. Such “sparse matrix” methods are used because the 
explicit inverse would be intractable or far more memory and 
computationally intensive than the sparse solution. To make 
conjugate gradient methods work quickly, a preconditioner 
must be selected that approximates the inverse solution to the 
true problem. While a solution can be obtained if no precon- 
ditioner is used (the identity preconditioner), the conjugate 
gradient method is greatly enhanced by the choice of precon- 
ditioner. The conjugate gradient algorithm is well-known in 
the literature and so we just review its application to the 
inverse of this problem. 

The problem we solve is (A /7 W 7/ WA+ e I)}(;=A /7 W // WI for 
the CARS process and $ H W H WS+€l)x=S H W H m for the 
CSRS process. These are total-least- squares solutions with 
Tikhonov regularization and weighting of the anti-Stokes 
(Stokes) spectrum. While the total -least-squares solution in 
general is not simply computable using simple operators such 
as the discrete Fourier transform, we can find an approximate 
inverse that is. We precompute y^=A H W H Wl and call the 
operator B=A H W H WA+€l so that this system becomes 
By=yJ for CARS. (For CSRS yC=S H W H Wl and 
B=S jf7 W 77 WS+€l). We assume that the anti-Stokes (Stokes) 
and Raman spectra are sampled to discrete frequencies. The 
matrix W is a weighting operator implemented by multiply- 
ing each anti-Stokes (Stokes) frequency by its corresponding 
weight. To speed up the solution, rapid ways of performing 
the operators 


3° 

A h , A, $> h , and S must be found. Fortunately, these can be 
computed using the Fast Fourier transform. To do this, we 
assume that the vector E is a sampled version of the incoming 
electric field E z (co) in the frequency domain. Likewise, we 
5 assume that R is a frequency sampled version of the reference 
pulse R(oo). We use the operator DFT to denote the discrete 
Fourier transform operation and DFT -1 to indicate its inverse, 
which is typically implemented using the Fast Fourier Trans- 
form algorithm. 

To compute A, we find the power spectrum of the instan- 
taneous intensity of the signal I £ =DFT[ I DFT - 1 [E] 1 2 ] . The 
operator Ax=(R*)DFT[DFT _1 [E]DFT -1 [1 E %\] . The discrete 
Fourier transform here is used as a rapid way to implement the 
15 convolution and product operations of Eq 6. Likewise, S %= 
(R*)DFT[DFT _1 [E]DFT _1 [I £ )(] !5! J (notice the complex con- 
jugate operator). Written this way, the Hermitian adjoints are 
easy to find. A^H^DFTLDFT-^EPDFT-^RI]] and 
S^H^DFTfDFT - 1 [E] *DFT“ 1 [RI] ] * . To save additional 
20 computation, the real-to-complex DFT can be used. 

Based on the forward operators for A and S, we can 
approximate an inverse for these operators that can be used as 
a preconditioner. These are based on a least- squares solution 
25 to each step of the forward operator, rather than computing 
the least-squares solution to the total operator. The approxi- 
mate inverse used for A and S are: 


30 


35 


/W ' |/ £ P + /|o i 


Rl 


\r\ 2 + i 2 0 


\DFr l [E]* 


| DFT 1 [E] I 2 + El 


h 

\Ie\ 2 + I 2 E0 


dft\ 


.r ri 1 . *r 

DFr 1 — 7 DFr l [E ] 

[\R\ 2 + I 2 ol 


I DFT ^ 1 [E ]\ 2 + El 


(9) 


(10) 


These are approximate inverses that have been regularized 
by three constants: I 0 2 , E 0 2 , and I^ 2 . These constants appear 
40 in divisors to prevent a division by a small magnitude number, 
which would produce instability in the solution. The con- 
stants are usually chosen based on the maximum value of the 
magnitude of the corresponding quantity. The constant I 0 2 can 
be chosen based on the maximum magnitude squared of a 
45 frequency component in the measured signal vector I, typi- 
cally some small fraction of it (e.g. 0.01). The constant E 0 2 is 
chosen by the maximum magnitude squared of a temporal 
sample of the sample field DFT _1 [E] (because the corre- 
sponding operation occurs in the time domain rather than the 
50 frequency domain). Finally, the constant \ EQ 2 is chosen by the 
maximum magnitude squared of any frequency sample of \ E . 

A summary of how to apply the preconditioned conjugate 
gradient algorithm to the inverse problem is given in FIG. 16. 
55 The method works by taking an initial random guess for the 

Raman spectrum that can be a random vector, a vector of 
constant values, or perhaps a prcondictioned conjugate gra- 
dient. From this vector the residual is calculated and reduces 
each iteration, with the two conjugate search directions 
60 updated. For this algorithm to work with CSRS rather than 
CARS, substitute S, S H , S approx ~ l for A, A H , and k appm ~ l in 
FIG. 16. 

Sometimes to improve accuracy, a conjugate gradient 
“restart” will be needed. To do this, the best solution from the 
65 last series of conjugate gradient iterations can be used as the 
starting point for the next series rather than a random guess. 
The best solution is found when the solution y} l) appears to 



US 7,586,618 B2 


31 

converge and stabilize near the solution, and then diverges or 
oscillates around the solution. Usually only one restart is 
needed. 

Stimulation of a Band of Raman Frequencies 

As an example of a pulse combination that can stimulate a 5 
band of Raman frequencies [2], we present a simulation of 
this technique using experimentally realistic values. We 
assume that the laser source is an ultrabroadband Ti-sapphire 
laser producing transform-limited pulses of uniform power 
spectral density between 700- 1 000 mil. The bandwidth from 10 
800-1000 nm is reserved for stimulating CARS, while the 
bandwidth from 700-800 nm is used as the reference pulse. A 
hypothetical molecule was used as the sample, with Lorent- 
zian resonances at 800 cm -1 , 900 cm -1 , 1000 cm -1 , and 1100 
cm -1 . These are frequencies in the Raman fingerprint region 15 
which is useful for practical molecular identification. A pulse 
with Q z =700 cm -1 , Q /7 =1300 cm -1 , and T=5 ps was used in 
Eq. 33 as the hypothetical stimulation pulse. The and 
are chosen to span the range of desired measured Raman 
frequencies. In practice, there should be extra measured 20 
bandwidth around the desired range (typically 20-25%) to 
ensure that the desired frequencies are sufficiently sampled. 
The time T is chosen to be on the order of the lifetime of the 
resonances to maximize resolution and signal. 

The results of the simulation are presented in FIG. 25 . Part 25 
(a) shows the magnitude of the illumination pulse, including 
the overlapped region where the Raman excitation occurs. 
Part (b) shows the spectrum of the illumination pulse. Part (c) 
shows the spectrum of the generated anti-Stokes light, simu- 
lated using Eqs. 1 and 2. Only the portion of the spectrum 30 
shorter than 800 nm is available for the inverse, because 
anti-Stokes in the same frequency band as the illumination is 
inseparable. Part (d) is the beat frequency spectrum, or the 
Fourier transform of the intensity of the pulse. This is the 
possible Raman frequencies the pulse can stimulate. Part (e) 35 
is the magnitude of the Raman susceptibility, and part (f) is 
the recovered Raman susceptibility using the regularized 
least-squares solution. 

A second simulation demonstrates the ability to recover a 
small portion of the Raman spectrum in detail. In this simu- 40 
lation, we desire to simulate the measurement of the relative 
amounts of DNA (deoxyribonucleic acid) which has a reso- 
nance at 1094 cm -1 , and RNA (ribonucleic acid) which has a 
resonance at 1 1 01 cm - 1 . A hypothetical Raman spectrum was 
created which has Lorentzian resonances at both frequencies . 45 
A pulse is designed such that Q z =1070 cm -1 , Q^ = l 130 cm -1 , 
and T=5 ps. The simulation of the anti-Stokes signal and the 
reconstruction of the Raman spectrum are shown in FIG. 26. 
Because of the relatively short pulse interval, there is a notice- 
able resolution loss to the spectrum, but the resonance lines 50 
are distinct. The short dephasing time in liquids, especially 
water is the practical limit on the instrument resolution in 
vivo. 
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What is claimed is: 

1. A method of examining a sample, comprising: 

exposing of the sample to a pump pulse of electromagnetic 

25 radiation for a first period of time, defined as the pump 
pulse duration, 

exposing of the sample to a stimulant pulse of electromag- 
netic radiation for a second period of time that is shorter 
than, or equal in duration to, the pump pulse duration, 

30 thereby producing a signal pulse of electromagnetic 
radiation for a third period of time which extends beyond 
the stimulant pulse duration, 

interfering the signal pulse with a reference pulse of elec- 
tromagnetic radiation, and 

35 demodulating a portion of the signal pulse that is not con- 
terminous with the stimulant pulse duration. 

2. The method of claim 1, wherein the sample is selected 
from the group consisting of a tissue sample, a single cell, and 
a patient. 

40 3 . The method of claim 1 , further comprising producing the 

reference pulse by exposing a reference material to a fourth 
pulse of electromagnetic radiation. 

4. The method of claim 1, wherein the first period of time 
is at least twice as long as the second period of time. 

45 5. The method of claim 1, wherein the electromagnetic 

radiation of the pump pulse and the stimulant pulse comprise 
electromagnetic radiation have a frequency in the range of 
infra-red to ultraviolet light. 

6 . The method of claim 1, further comprising producing the 

50 pump pulse and the stimulant pulse from a single laser. 

7. The method of claim 1, further comprising producing 
digital data from a portion of the signal pulse not produced 
during the exposing of the sample to the stimulant pulse. 

8. The method of claim 1, further comprising producing an 

55 image from the digital data. 

9. A method of producing an image, comprising collecting 
a set of data, wherein the data comprises a plurality of digital 
data produced by the method of claim 7. 

10. The method of claim 9, wherein the image is an OCT 

60 image. 

11. The method of claim 1, wherein the pump pulse is 
chirped. 

12. The method of claim 6, wherein the laser has thousands 
of cm -1 of bandwidth. 
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